Formation of orbital-selective electron states in LaTiOa/SrTiOa superlattices 
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The interface electronic structure of correlated LaTiOs/SrTiOs superlattices is investigated by 
means of the charge self-consistent combination of the local density approximation (LDA) to density 
functional theory (DFT) with dynamical mean-field theory (DMFT). Utilizing a pseudopotential 
technique together with a continuous-time quantum Monte-Carlo approach, the resulting complex 
multiorbital electronic states are addressed in a coherent fashion beyond static mean-field. General 
structural relaxations are taken into account on the LDA level and cooperate with the driving 
forces from strong electronic correlations. This alliance leads to an Ti(3da;j/) dominated low-energy 
quasiparticle peak and a lower Hubbard band in line with photoemission studies. Furthermore 
correlation effects close to the band-insulating bulk SrTiOs limit as well as the Mott-insulating bulk 
LaTiOs limit are studied via realistic single-layer embeddings. 

PACS numbers: 73.20.-r, 71.27.-|-a, 71.15.-m 



Research on layered hcterostructures composed of 
various different transition-metal oxide (TMO) com- 
pounds emerge as a major new field in condensed mat- 
ter physics pi Especially the intriguing physics of a re- 
sulting two-dimensional (2D) electron gas from interfac- 
ing bulk-insulating TMOs may open the door for new 
tailored hybrid materials with specific transport, mag- 
netic and/or superconducting properties. Albeit there 
are various such different layered TMO combinations re- 
alized, heterostructures build up by interlacing band- 
insulating with Mott-insulating compounds are partic- 
ularly appealing. As they raise questions about the elec- 
tronic states resulting from conceptually quite different 
limits, this structural setup challenges the existing mod- 
ern first-principles approaches to electronic structure. 
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Coherent superlattices (SLs) of the strongly correlated 
LaTiOa (LTO) Mott insulator with the SrTiOa (STO) 
band insulator belong to the most prominent examples 
of these structured materialsi^ Not only displays the 
LTO/STO interface metallicity, also magnetic correla- 
tions are vital due to the antiferromagnetic ordering of 
LTO^ below rN=146K. Moreover 2D superconductivity 
was revealvedi^ Aside from possible polarization mech- 
anisms^ and suppressing Mott correlations, the nomi- 
nal Ti^+ valence in LTO allows for doping the band- 
insulating STO side, giving rise to metallic transport. 
Yet in reality strong Coulomb interactions among the 
electrons complicate this simplistic picture. Photoemis- 
sion experiments^ indeed reveal strong correlation signa- 
tures, i.e. lower Hubbard band and quasiparticle (QP) 
peak, in the valence spectrum. The interface conduc- 
tivityi as well as the optical response^ is identified to 
dependent on electronic correlations. Numerous theo- 
retical studies on the interface electronic structure exist. 
There are inital tailored Hubbard-model considerations 
from unrestricted Hartree-Fock;^ several first-principles 
investigation based on Kohn-Sham density functional 
theory (DFT)^-i^ DFT-^U studiesHOi ^g well as many- 
body approaches based on the Lanczos-mcthodp^ slave 
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FIG. 1. (color online) Left: (4,4) x2 LaTiOs/SrTiOa 
ideal tetragonal unit cell with La (big green/grey), Sr (big 
blue/dark), Ti (grey/light grey) and O (small red/dark) ions 
as well as the charge density from t29-like LDA bands below 
ep, holding 8e~. Right: corresponding LDA bandstructure 
and DOS for ideal-tetragonal and orthorhombic-relaxed case. 



bosonSii^ and dynamical mean-field theory (DMFT) i^^d^ 
Contrary to former studies, this work treats the ef- 
fective single-particle character of the materials chem- 
istry on an equal footing with many-body effects from 
a local perspective. We performed charge self-consistent 
DFT-HDMFTi^i^o computations for selected LTO/STO 
SLs that allow for detailed examinations of the subtle 
interplay between realistic interface effects and multi- 
orbital electronic correlations at room temperature. It 
is revealed that structural relaxations and electronic cor- 
relations ally in driving an enlarged Ti{3dxy) orbital po- 
larization in real space and in the low-energy spectrum. 
Fostered by enlarged lateral coherency effects, a promi- 
nent renormalized d^y QP peak resides close to the Fermi 
level. Moreover investigated single-layer LTO/STO ar- 
chitectures remain metallic troughout the SLs. 

The theoretical approach (see Refs. [2l| andl22l for de- 
tails) builds up on the combination of a mixed-basis pseu- 
dopotential framework^^ with a hybridization-expansion 
continuous-time quantum Monte Carlo solver—^— for the 




FIG. 2. (color online) (a) Inequivalent Ti ions for (4,4) x 2. (b) Local Til-5 LDA DOS for unrelaxed (grey background (bg)) 
and structurally relaxed (white bg) cases, (c) Local DFT+DMFT spectral function for Til-5 with same color coding. Though 



dx 



are different by symmetry with relaxations, Z and n differ only marginally, hence averaged values are shown. 



DMFT impurity problem. To include important struc- 
tural relaxation&iir— in a general scope, allowing for 
layer-distance variation and tilting of the TiOg octahe- 
dra, we constructed SLs in an (n,m) x 2 setup, where n, m 
denote the numbers of LaO, SrO layers in the unit cell 
with two inplane Ti ions, respectively (see Fig. [T]). The 
two Ti ions are assumed equivalent by symmetry in each 
layer, so possible lateral orderingsi^ are suppressed. For 
all discussed SLs the lateral lattice constant was set to 
the STO value a=3.905A and a ratio c/a=0.99 was iden- 
tified reasonable. 

Figure [T] shows the LDA bandstructure for a (4,4) x 2 
superlattice (80-atom unit cell) along with the density of 
states (DOS). The occupied bands just below the Fermi 
level Sf with dominant Ti(i2g) character accommodate 
the eight additional electrons from the La'^^ ions in the 
unit cell. Structural relaxations, relevant also in the STO 
part, enlarge the gap between the 0(2p)-derived bands 
deep in energy and the latter t2g bands as well as in- 
crease the pseudogap between t2g and gg-like states high 
in energy. A DOS maximum right at sp exists for both 
structural cases, but there is a gain of 36 meV/atom in 
the LDA total energy upon relaxation. The real-space 
distribution of the occupied i2g valence charge density in 
Fig. [T] elucidates the t2g doping in the STO part. This 
charge transfer is even strengthened in the relaxed or- 
thorhombic structure. 

The (4,4) X 2 unit cell contains five symmetry inequiv- 
alent Ti ions, denoted here Til-5, which local t2g DOS 
from projected local orbitals^ is displayed in Fig. [2]3. 
From the mid LTO part (Til) to the mid STO part 
(Ti5) the t2g filling is decreasing. The effective band- 
width shrinks from ^4 eV down to ~2.5 eV, with mi- 
nor smaller size in the relaxed structure. Though over- 
all rather balanced, the dxy occupation is somewhat in- 
creased by structural relaxation. In the latter case, the 
total number of t2g is higher, i.e. the doping of these 



states is more efficient. 

To capture the effect of many-body correlations an ef- 
fective three-orbital Hubbard Hamiltonian T-L with fully 
rotational invariant interaction terms, i.e. 
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is applied at each individual Ti site i. Its parametrized 
by the adequate^LS^ Coulomb integral f/=5eV and 
the Hund's exchange J=0.7cV with U'=U—2J and 
U"=U~3J. Including symmetry, this leads to five in- 
equivalent single-site impurity problems embedded in the 
full charge self-consistent DFT-I-DMFT calculations for 
the (4,4) X 2 unit cell. Charge self-consistency is a vi- 
tal methodological ingredient because of the subtle elec- 
tron transfers (cf Fig. [1]). For the projected-local-orbital 
construction of the correlated subspace a number of 80 
Kohn-Sham bands starting from the bottom of the ^23- 
likc manifold was used. The double-counting correction 
applied to each Ti impurity self-energy Ej amounts to 
an site-avcragcd fully-localized^^i term. All DFT-hDMFT 
computations were performed at r=290K^ 

In Fig. [5J; the resulting local t2g spectral functions are 
plotted. Compared to LDA the total filling increases 
once more with correlations. Close to the interface es- 
pecially the dxy orbital character gains further occupa- 



tion with correlations. Moreover the d, 
e-p is now pronounced compared to dx 



xy 
z,yz, 



weight near 
resulting in 

a dominant dxy QP peak below the Fermi level. Coher- 
ent transport is thus dxy dominated. From the incoherent 
high-energy part, the correlation strength is larger in the 
structurally relaxed orthorhombic SLs. Right at the in- 
terface (Ti3 ion) the dxy local spectral part is susceptible 
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FIG. 3. (color online) Total t29-like DFT+DMFT spec- 
tral function compared to the LDA DOS for (4,4) x 2. Top: 
tetragonal unrelaxed (grey bg), bottom: relaxed orthorhom- 
bic (white bg). Right part shows a blow up close to ep. 



to (pscudo)gapping. The site- and orbital-resolved QP 
weight Zim — {^~dY.im/duj)~^ varies significantly accross 
the interface, revealing somewhat stronger mass renor- 
malization for d^^y. The complete DFT-^DMFT t2g-like 
spectral function in Bloch space carries the dominance 
of the dxy close to the Fermi level (see Fig. [3]) . While in 
the unrelaxed tetragonal case a maximum at ep remains 
as in LDA, a minimum occurs when including structural 
relaxations. Increased spectral-weight transfer towards 
the lower Hubbard band takes place in the latter. This 
Hubbard peak shifts closer to —1 eV with relaxations, in 
good agreement with photoemissioni^ 

In order to obtain deeper insight in the relevance of 
correlation effects, lets turn now to the limiting case of a 
single LaO(SrO) layer within an STO(LTO) host. Within 
our superlatticc approach two different unit cells were 
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FIG. 4. (color online) Left: total fos-like DFT-f DMFT spec- 
tral function compared to the LDA DOS for (1,9) x2 (top) 
and (9,1) x2. Right: Ti(i29) occupations with distance from 
the LTO/STO layer within LDA (left part) and DFT-I-DMFT 
(right part). Circles mark results for the 4-layer host, trian- 
gles for the 9-layer host. 



chosen, respectively, incorporating 4 and 9 host layers. 
The larger structure amounts to a 100-atom unit cell. 
While the 4-layer-host case has 5 Ti layers inbetween 
the single layers and 3 inequivalent Ti ions, the 9-layer- 
host structure has 10 Ti layers and 5 inequivalent ones, 
denoted Til-5 with increasing distance from the threaded 
single layer. Hence the smaller(larger) structure has an 
odd(evcn) number of Ti layers. In the following only the 
structurally relaxed orthorhombic cases are discussed. 

Figure |4] displays the spectral comparison between 
LDA and charge self-consistent DFT-f-DMFT applied to 
the (l,9)x2 as well as the (9,l)x2 structure. The cor- 
related spectral function of (1,9) x 2 with one LaO layer 
in STO is metallic in agreement with optics;^ showing 
a smaller QP peak below ep and a more prominent 
one just above. While at low energy the LDA differ- 
ence between both structural types amounts mainly to a 
Fcrmi-level shift (1 vs. 9e- below ep), DFT+DMFT sig- 
nals the increased correlations for (9,1) x 2 via substantial 
spectral-weight transfer to Hubbard bands. No insulat- 
ing state, also not site-selective (see below), is obtained 
for (9,1) x 2, which thus could be viewed as a doped Mott- 
insulator. Note that also bulk Sri-ajLa^jTiOa is insulat- 
ing only above a;c^0.95i^ Experiments on confined STO 
in a larger GdTiOs host yet show the principle chance 
for reaching an insulating interface i^ 

Instructive are the local t2g-like occupations for Til- 
5. For the single LaO layer the Ti doping in STO is 
still not accomplished in (1,4) x2 far from the interface. 
But for (1,9) x 2 the Ti5 ion indeed shows zero t2g fill- 
ing (see also FigE]). DFT+DMFT leads here to some- 
what farther charge flow into STO and a stronger orbital- 
filling alignment within the t2g manifold away from the 
interface. However again the d^y polarization close to it 
is strengthened with correlations. The evaluated sheet 
carrier densities (in unit cm^^) n^[^^j=1.8xl0^* and 
n^^J+°'^'^'^=2.3xlO" agree well with the experimental 
value "shcot^"^^-'^'-'^^ from optics4 The increased charg- 
ing of the Til ion in DFT+DMFT for the single SrO 
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FIG. 5. (color online) Left: LDA charge density from occu- 
pied t2g-like bands for a Ti column along the c axis of (1,9) x2. 
Right: same for (9,1) x2. Middle: LDA t2g-like DOS (left) 
and DFT+DMFT spectrum (right) for Til, Ti5 (see text). 




FIG. 6. (color online) Formation energy and convex hull for 
the embedded single-layer SLs. 



layer architecture (again favoring d^.y) is evident. In gen- 
eral, whereas LDA quickly saturates here to the nominal 
n=l t2g occupation, many-body effects result in a bal- 
ancing of the strong LDA occupation differences with 
distance from the interface. Even for Ti4,Ti5 the hole 
doping is vital. A subtle d^y/dyz filling crossover occurs 
near Ti2, which marks the competition between bulk- 
LTO Mott-insulating behavior and LTO/STO interface 
physics. While in the former case indeed the dyz orbital 
has dominant contribution to the correlated crystal-field 
ground statcj^^ the driving force for dxy polarization is 
stronger at the interface. 

Besides illustrating the real-space variation of the low- 
energy LDA valence charge density for (9,1) x 2 and 
(1,9) x2, in Fig.[5]the local t2g-hke spectral properties are 
visualized for the Ti ions closest (Til) and farthest from 
the single LaO(SrO) layer. The QP structure for Til 
in (1,9) x2 is more subtle than for the interface ions in 
(4,4) X 2, since a strong orbital-selectivity towards dxy has 
not yet manifested in the single-layer set up. Still a shift- 
ing of the large d^z.yz DOS below ep towards the unoccu- 
pied part is visible. The i2g-hke orbital behave more bal- 
anced in (1,9) X 2 with significant Hubbard-band weight 
for Ti5. However a standard doped- Mott-insulating pic- 
turing of the Ti5 multiorbital spectral function seems not 
applicable, i.e. the QP weight remains sizeable. 

Finally we want to comment on the energetics of 
the relaxed embedded single-layer SLs within LDA and 
DFT+DMFT3^ To this we compute the formation en- 
ergy, defined here as 



Ch&~n/{n+m). Be aware of the nontrivial character, 
since in the correlated case not only i?ff'™ includes 
many-body corrections but also E^q grpQ . As the pro- 
jected local orbitals are here derived from f2g-like bands, 
however that correction vanishes for bulk-STO having 
those unoccupied. Yet bulk-LTO is a Mott insulator in 
DFT-f DMFT^ which is verified within the charge self- 
consistent scheme for the here given many-body Hamilto- 
nian ^, chosen Coulomb parameters and double count- 
ing. Within our double-counting scheme the bulk-LTO 
correlated Mott state enters eq. ^ with a much lower to- 
tal energy. Figure [6] shows the variation of the formation 
energy based on the four single-layer (n,m) x2 structures. 
Already in LDA an asymmetry of the basic convex hull 
towards the LTO side exists, which is enforced with cor- 
relations. Given the delicacy of E^°'^'^, the values within 
both theoretical schemes are rather similar, which espe- 
cially speaks for the novel DFT-I-DMFT approach to such 
subtle energetic quantities. The observed trend of gen- 
erally reduced 5"^°™ (though evaluated at finite T) may 
correct frequent LDA overestimations. 

In summary, the correlated electronic structure of re- 
alistic LTO/STO SLs has been studied witin charge 
self-consistent DFT-I-DMFT beyond static mean-field. 
Many-body effects and structural relaxations are revealed 
to ally in driving orbital-selective behavior towards dom- 
inant Ti{3dxy) filling and QP behavior. In addition, sig- 
nificant spectral-weight intensity just above the Fermi 
level renders intricate susceptibility to applied field pos- 
sible. Single-layer architectures of these systems result in 
global, but still intriguing, metallic behavior for the stud- 
ied embeddings. In general, the recent advances put the 
DFT+DMFT formalism in position for challenging ma- 
terials investigations, including engineered systems. A 
further methodological step will be the sound computa- 
tion of corrclation-infiuenced structural relaxations. Fu- 
ture DFT+DMFT work on TMO hetcrostructurcs such 
as the study of possible ordering instabilities, of different 
interface geometries and of spin-orbit effects is envisaged. 
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